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GRAPHIC DESIGN OF ALLUVIAL CHANNELS 


Ning Chien,’ A.M., ASCE 


SYNOPSIS 


The determination of channel depth and slope to carry a given unit dis- 
charge and sediment load is given in graphical form on the basis of the 
Einstein bed-load function. Applications of these curves are made to evaluate 
the effects on river regimen by the diversion of flows, construction of dams, 
contraction of channels, elimination of river forks, and cut-off of bends. The 
deficiencies of knowledge in describing the canal flows and in determining the 
sediment supply to the canal systems still prevent the use of the method in 
designing irrigation canals. 


INTRODUCTION 


The design of alluvial channels to carry a given flow is one of the impor- 
tant problems in the field of hydraulic engineering. As early as in the 
eighteenth century, Chezy developed a formula which relates the depth and 
slope of the channel to the velocity of flow. As time passed, engineers soon 
realized that the problem is complicated by the fact that alluvial streams 
carry not only the water flow, but also sediment of the same type as that 
which composes the bed. Erosion or deposition will occur, if no proper pro- 
visions are made for the sediment flow, and thereby change the character- 
istics of the channel completely. This paper attempts to show how one can 
utilize the results of sediment research to develop a more rational approach 
to the problem of alluvial channel design. 


The Problem 


The problem of alluvial channel design can be stated as follows: Given the 
bed material, what should be the width, depth, and slope of the channel to ac- 
commodate the given discharge and sediment inflow? The sediment inflow 
and the discharge can be regarded as conditions imposed on the channel by 
the watershed and are independent of the channel itself. For a channel in 
equilibrium, we have two conditions to be satisfied and three variables. 

This leaves one of the three variables remaining undetermined, according to 
the present state of knowledge. 

It is to be noted that the sediment inflow refers essentially to the bed- 
material load (EINSTEIN and JOHNSON, 1950). The wash load, in spite of its 
often being the major part of the sediment load and the deciding factor in the 
silting of reservoirs and lakes, plays only a minor role in the problem of 
channel stability. Indirectly, the presence of large quantities of wash load 
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may be a determining factor in the use of permeable structures for construc- 
tion works and regulating measures. 


The Design Curves 


The design of alluvial channels is greatly facilitated by the use of a graph- 
ical presentation first suggested by RAYNAUD (1953). The theories de- 
veloped by EINSTEIN (1950) in describing the alluvial channel flow and in re- 
lating the sediment transport rate to the flow conditions are taken as the 
basis for the construction of the design curves. For illustrative purposes, 
the following simplifications are introduced such that the curves are of more 
general use: 

a) bank friction is neglected; and 

b) bed material is uniform in size. 

It is understood that in the actual design where there are not more than a few 
particular channels to be considered, one can always take all the local condi- 
tions into consideration and thereby construct more reliable curves. 

Figs. 1-3 give the design curves for the determination of the slope and 
depth with a given discharge and sediment inflow (including both the sus- 
pended-load and bed-load) for channels of fine sand (D = 0.25 mm), of coarse 
sand (D = 2.5 mm), and of gravel (D = 25 mm), respectively. The discharge 
and sediment inflow are expressed in terms of quantities per unit width of the 
channel, leaving the variable width open. These curves are applicable only 
up to the bank full stage. 

The absence of the width from the design curves is not as detrimental as 
it first appears. Plans on river improvements usually include a certain 
amount of bank stabilization work, and the width of the channel can be chosen 
freely within limits. In nature, river channels may also be confined in 
canyons, or the banks may be consolidated by the growth of vegetation. On 
the other hand, LACEY (1939) claimed that the width, p, of a natural channel 
is related to the discharge, Q, by the following formula: 


p = 2.67 Q'/2 (1) 


LEOPOLD AND MADDOCK (1953), after analyzing many streams in the Great 
Plains and the Southwest of the United States, have come to the conclusion 
that (a) the relationship between channel width and discharge as given in eq. 
(1) is valid only for the comparison of various cross sections at correspond- 
ing discharges; (b) the constant of proportionality in eq. (1) is actually, among 
others, a function of the bed material size; and (c) the channel width varies 
only with the 0.26 power of the discharge when considering the variation in a 
particular cross section of a river. 

An examination of Figs. 1-3 indicates that for a given slope (or depth), the 
sediment transport rate increases substantially with a slight increase in 
depth (or slope). This implies that so far as the design of alluvial channels 
is concerned, considerable tolerance can be allowed in the determination of 
the sediment transport rate, a fact which should be welcomed by engineers 
who may have doubts in their minds concerning the reliability of existing 
sediment transport theories. Figs. 1-3 also show that an increase in dis- 
charge results in a much larger increase in its capacity to transport sedi- 
ment. This has considerable bearing on the planning of river works, and will 
be made clear in the examples cited in the next section. 
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Planning of River Works 


In planning river works, the problem is often not how to design a new 
channel to carry the discharge and sediment load, but rather how to regulate 
and improve the existing channels such that they will function properly. Any 
alteration of the channel naturally upsets the existing regimen of the river; 
and it is very important to know how the channel adjusts its width, depth, and 
slope to achieve the new equilibrium. A comparison between the existing 
regimen and the future regimen will then reveal whether or not the plan 
serves its purposes. In this respect the use of design curves like Figs. 1-3 
is of great help. In the discussion to follow, hypothetical examples will be 
cited to illustrate how different engineering measures affect the character- 
istics of streams. 


1) Diversion of Flow 

Flow is diverted from streams either to relieve downstream flooding, or 
for irrigation and other domestic and industrial uses. A division of flow 
naturally induces certain changes downstream, and it is important to ascer- 
tain whether the changes will lead to adverse consequences. For instance, 
for some time part of the Mississippi River flow has been diverted to the sea 
through the Old River-Atchafalaya channel 320 miles upstream from its 
mouth, The latter channel is enlarging at an alarmingly rapid rate, and cur- 
rently takes 23% of the Mississippi River flow. A control structure has been 
planned by the corps of Engineers to be installed in Old River. Unless this 
is done during the next 20 years, most of the flows will be diverted away 
from the lower Mississippi, leaving the port of New Orleans more or less 
isolated. A crisis of this kind arises essentially from an unfavorable divi- 
sion of the sediment load between the two channels, causing aggradation in 
one channel and degradation in another. The amount of sediment drawn by 
the branch channel depends very much upon the local conditions, the con- 
trolling factors of which have been clearly pointed out by THOMAS (1953). 
One must either control the division of flow according to the existing condi- 
tions, or change the flow patterns at the point of diversion by engineering 
measures such that the resulting load division will cause neither erosion nor 
deposition downstream. In either case it is the question of what should be the 
division of the sediment supply for a given division of flow that will result in 
equilibrium flow at existing channel slopes. This can best be answered by 
plotting the sediment rating curves at constant slopes, a few of which are 
given in Fig. 4 for different channels. The smaller the slope of the curves is, 
the larger is the percentage of sediment load that must be diverted by the 
branch channel without upsetting the equilibrium of the main channel down- 
stream. A few numerical examples are also drawn from Figs. 1-3 and given 
in Table I. The subscripts 1 and 2 refer to the channels upstream and down- 
stream of the division point, respectively. In case I no change in width is 
assumed to take place as the result of flow diversion, and in case II the width 
is assumed to change proportional to the 0.26 power of the discharge. 

An examination of Fig. 4 and Table I reveals the following: 

a) As the slope of the sediment rating curves is smaller than 45 degrees, 
the downstream channel will not aggrade only when a larger percentage 
of sediment load is diverted in comparison with the percentage of flow 
diversion. 

b) For the conditions given in Table I, a diversion of 50% and 25% of the 
flow calls for a diversion of 65-90% and 35-75% of the sediment load, 
respectively, if the channel width remains constant. These figures for 
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TABLE I 
EFFECT @ FLOW DIVERSIONS 


cfs/ft, tons/day 
| f/f 


___ Case I, B=const. 


. % is the percentage of the original flow which enters the downstream channel, 
Ql- S/%) is the percentage of flow which is diverted, 


** Qra/Qy, ie the percentage of the original bed-material load that is allowed to 
enter the downstream channel without changing the channel slope. 


the required diversion of sediment load are reduced to 60-85% and 30- 
1 55% if the channel width changes also with the discharge to the 0.26 
power. 

c) The percentage of sediment load that must be diverted to satisfy the 
original equilibrium conditions of the channel (i) decreases with the 
slope and the discharge for the same bed material size; (ii) increases 
with the bed material size for small streams which carry little sedi- 
ment load; and (iii) decreases with the bed material size for large 
streams with heavy sediment load. 


2) Construction of Dams 
The changes in river regimen caused by the construction of dams across 


the river have received considerable attention in recent years. The 
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mm, % (\cfs/ft. | day/ft | efs/ft. tons/day tn 
| 
| | b&b (11381 6 7 | 183 
| £00025 195 | | 1963 
| 125s 7 ; 100 19. 
75 | 100 | Bins | jms | 
| ,00025 195 100 | (7.5 | 
00050 | 1400 _|$0,.0 800 | 
300025 | _ 2200 000 | 1300 | 
| 265 | 100 | 2 ae, | 
| | 250 30 32.3 | 150 1590’ 
| 4 i j 
| 250 167.5 63.5 | 200 | 230 | 62,8 
+0015 7h0 | 500 62.7 
2000 1270 | 6365 1350 62.6 | 
25 50 i .003 13| 50 - | -| © 1 5.2 | 
005 | 1350 | 305 | 22.6 | 
2005 1350 760 56,3 850 58.5 
010 00 B00 66 500 2 


completion of a number of large dams, including the Hoover and the Elephant 
Butte Dam, also offers a rare opportunity for close field observations. Asa 
result, it has become common knowledge that damming up the flow will in- 
duce aggradation upstream and degradation downstream, two of the immedi- 
ate effects which can also be derived easily from Figs. 1-3. However, it 
becomes less obvious when comes the question of what will be the changes in 
the river regimen below the junction of a sizable tributary downstream from 
the dam. The elimination of peak flows by the regulation of the reservoir 
greatly reduces the capacity of the channel to transport sediment. In fact, if 
the total annual discharge is released from the reservoir at a constant rate, 
the downstream channel may carry less than one-tenth of the sediment load 
that the same total annual discharge carries if it flows according to its 
natural pattern. Or, in other words, the building of a dam upstream reduces 
both the sediment capacity and the sediment supply in the channel below the 
junction point, the net effect of which can only be determined by examining 
the design curves. 

Consider, for example, the junction of the two channels as shown in Fig. 
5. Will the channel below the junction degrade or aggrade if the discharge 
from the main channel is cut to 20,000 cfs by the regulation of a reservoir? 
Table Il(a) gives the characteristics of these three channels before the 
construction of the dam in channel (1). It can be seen that the slope of the 
channel below the junction is lower than that of either of the confluents, a 
fact substantiated by many field observations (MACKIN, 1948). The slope 
of the main channel below the confluence with a sizable tributary will be 
steeper than that above the junction only when (a) the tributary brings in 
sediment substantially coarser than that carried by the main stream, such 
as the increase in Missouri River slope below the mouth of the Platte River; 
or (b) the tributary carries a disproportionately large sediment load with 
respect to its discharge, such as the steepening of the Rio Grande River 
below the confluence of Rio Puerco River. Table I(b) indicates the change 
in channel characteristics after the completion of a dam in channel (1). 
Two cases are studied. In case I the degradation below the dam has not 
yet been extended to the reach immediately above the junction such that 
the slope of the channel in that reach remains unchanged. In the second 
case, the degradation below the dam is completed and channel (1) be- 
comes essentially non-alluvial in nature and carries practically no 
sediment. In both cases the channel below the junction is substantially 
accreted. This is true also for channels other than these cited in the 
example above. It appears that the effect of reduction in peak discharge 
by the operation of the reservoir by far outweighs the effect of decrease 
in sediment supply, and one may regard the aggradation of the main 
stem downstream of the dam and below the junctions with tributaries as 
a rule common to many alluvial streams. 


3) Contraction of Channel 

One of the most common measures in increasing the navigable depth 
of a river is to confine the channel between jetties or training walls for 
limited or extensive stretches. A reduction in width means a propor- 
tionate increase in unit discharge and unit sediment supply, the effect of 
which on the channel characteristics can be derived from the design 
curves. Table I shows the changes in slope and water depth if the 
channel is contracted by 50%. Generally speaking, this results in an in- 
crease in depth of 45-80% and a decline in slope of 15-30%. STRAUB 
(1934) has developed a formula expressing the depth in the contracted 


611-5 


€ 


| 
- 
|secoo 


Ts w 
| | edors 


uyeu Uy weg JO ey} edots uy eBueyo (T) 


( UTey weer qedg ey) UT weg P 


(q)II 


tpeazesgo St (C€S6T) pue pTodoe] Aq peuteygo 
~ (2) 


seat [90004 stri-loroz zz | att | | coo'set | | 

| | 6g000 = | 02000° ot os 00S ¢)000%S o00*sz 

| 9*9T | 99000 = jo*gt | oot | | Tt 


611-6 


| 
lg 
| 
Lt | 
| | abl alae 
a 
| 
| 
| PROG 
| 8 | i 
| | 
} 
| 
| slo; e 
j 
| 8 
: 
eat 1g 
8 
> 


TABLE III 


Effect on the Channel Characteristics 
by Contracting the Channel by 50% 


| 

002 12.5/ Ob | 6.75! 25.0 0.8 ,000157 | 10.5| -21.5 | 55.6 | 
25.0 2.5 | 8.80. 50.0 ; 5,8. 2000150 13.9 | 58.0 


| $040 | 1245 | 10,75 100.0 | 17.4 -27.5 | 61.8 | 
2.90, 10.0 0.9; 2000300 4.6) -25.0 60.0 


| 
| 100:0 80.0 | 13.30 200.0 160.0 .0001L0 | 21.0 | -30.0 | 57.9 


3.90 25.0 8.0 4000300 | 6,2 =25.0 | 59.0 
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6.50 100.0 250.0 | 2000285 10.6; -28.7 } 63.1 
8.95 200.0 1200.0) .000230/ 16.5| | 
| 13.30 8000.0 2000255 | 23.6) -36.2 | 77.5 
3225, 25,0 | 8.0) | | 7h.5 
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6.20! 100.0 . 250.0 0017 9.2 -15.0 48.3 
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| 
0.1 | 5.40 25.0 C,2 200073 


7.00 50.0 4.0 00072 12.0 71.5 
9.65) 100.0 30.0 .00077 -23.0| 59.6 
10! 65.0| 12,80; 200.0| 130.0) .00083 18.5 -17.0; 
200.0! 180.0 | 16,00! 400.0 | 360.0 .00075 26,0) -25.0 | 62.5 

50.0) 0.2) 9.90! 100.0 0020 | 18,0| -33.3! 61.8 | 
100.0 H 12.5 14.20, 200.0 25.0. 0022 23.1 62.7 
200.0! 140.0! 280.0} .0023 -23.3 76.2 

25.0 1.5| b.70} 50.0 3.0| | 8.1/ -30.0 72.3 
50.0 20 6.60 100.0 80.0 11.7 -31.6 77.3 
100.0 300.0 9.” 200.0 600.0 20047 14.9 -21.6 
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section in terms of the flow conditions in the uncontracted section as follows: 


(2) 4d, (4, 


+1 
d, = -a EF 1-@a@ (2) 


in which water depth at the contracted section, 

water depth at the un-contracted section, 

critical tractive force, 

slope in the un-contracted section, 

unit weight of the water, 

the amount of contraction 

(p, - P.)/P,, P, and p, being the widths of the channel at the 
normal and contracted sections, respectively. 


For most of the streams, the critical tractive force is small and can be 
ignored. Eq. (2) thus becomes 


d, = d,/(1 - (3) 
and 
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(da, - d,)/d, = 1/(1 - 1 


a equals 0.5 for 50% reduction of channel width; and eq. (4) gives an in- 
crease in depth of 56%, a value which checks well with the results given in 


Table I. 


4) Closure of River Forks 
When a river flows in a divided channel, the water depth in either of the 
channels is often insufficient for navigation at low water, and it is desirable 
to eliminate the forks. Consider, for instance, half of the flow goes down 
each arm as shown in Fig. 6: Table IV shows the channel characteristics of 
the river at the fork and in the undivided reach. The water depth in the forks 
is only 80% of the depth in the single course if the flow is equally divided at 
the fork. It is interesting to note that FRANZIUS (1936) gives the ratio of 
d,/d, from 0.777 to 0.80 by making several assumptions concerning the hy- 
draulics of the channel. 


TABLE IV 


Channel Characteristics of a River 
At the Pork and in the Undivided Reach 


Undivided Channel Chmnels at the fork 


5) Cut-Off 

When a stream flows in excessive meander loops, the water at high stages 
may take a short-cut across the narrow necks of overdeveloped bends. This 
natural process is often assisted by engineering measures so as to relieve 
the flood pressure upstream and, to a lesser degree, to increase the maneu- 
verability of barge fleets. A cut-off of the bends means a shortening of the 
river course, and consequently, a local steepening of the channel. This local 
increase in slope will eventually be wiped out by lowering of bed upstream 
and by a raise of river bottom downstream. If the channel is subject to slow 
accretion, the step will be absorbed in the final adjustment of the profile, 
reducing the aggradation of the upstream channel to an amount equal to the 
product of the length of the cut-off channel and the final equilibrium slope. 
The reliance on cut-offs in reducing flood stages upstream, therefore, is 
only effective when either the channel is slowly building up at the present, or 
the downstream aggradation causes no serious local damages. 

It has been demonstrated so far that valuable information on the inter- 
dependent effects of changes of the channel characteristics, of the sediment 
supply, and of the flow can be derived from proper use of the design curves. 
In addition to this, there are other complications which must be considered 
so as to obtain a general picture as complete as possible. Any change in 
conditions may lead to other changes, and thereby further complicate the 
final adjustment. For instance, lowering of the channel bed may increase 
the sediment inflow from tributaries resulting from the increase in slopes 
near their mouths (LANE, 1942). There are also conditions which are 
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peculiar to the particular stream under consideration. No two naturai 
streams are exactly alike. Each has had its own complex history, and 
possesses certain features and irregularities not to be found in other similar 
drainage basins. For instance, the existence of gravel lenses or outcrop of 
rock in the bed may allow the bed to be scoured only down to a certain extent, 
beyond which the flow is forced to widen the channel which may not be the 
case otherwise. An adjustment of the slope may be accomplished by the dif- 
ference in channel filling or scouring between the upstream and downstream 
part of the reach, or may be achieved by the change in meander pattern, de- 
pending very much on the erodibility of the bank material. Rubey (1952) has 
also pointed out that a change in regimen is liable to be counter-acted by an 
adjustment in cross-section in some cases, and by modification in slope in 
others, depending on the historical background and physiographical and geo- 
logical features of the watershed. 

The design curves also do not give any clue on the time factor involved in 
planning the river works. Besides the outcome of the regulating measures, 
it is just as important to know how long the river takes to reach its new 
equilibrium in response to the changes. A plan which will bring immediate 
relief, but in the long run will retreat to its old conditions, may still prove 
to be economically feasible if the time involved between the immediate change 
and the final outcome is sufficiently long. 


Design of Irrigation Canals 


The design curves as given in Figs. 1-3 cannot be used for irrigation 
canals design because of the lack of certain information. The difficulties lie 
partially in the description of the hydraulics and partially in the determina- 
tion of the sediment supply. 

Irrigation canals are usually of much smaller size than natural streams. 
The bank friction is a significant part of the total resistance. The amount of 
bank friction depends very much on the shape of the channel, thus making the 
construction of generalized design curves impossible. Lacking are also sys- 
tematic studies on the factors which govern the additional resistance caused 
by the existence of sand bars and other irregularities at the bed surface of 
narrow channels. It has been found for natural streams that the bar resist- 
ance is a function of bed-load transport (EINSTEIN and BARBAROSSA, 1952), 
and the functional relationship between the two has been incorporated in the 
preparation of Figs. 1-3. For more confined canal flows, this empirical re- 
lationship between the bed-load transport and the bar resistance appears to 
be somewhat different. There are also indications that other factors, such as 
the width of the channel, may also come into the play. This deficiency of 
knowledge calls for prototype measurements which must be analyzed in the 
light of modern concepts on friction in open channel flow (EINSTEIN and 
BANKS, 1951). 

The determination of sediment inflow in canals is rather complicated. In 
planning river works, the sediment supply to the reach under consideration 
can be calculated by applying the Einstein bed-load function to the reach im- 
mediately upstream. The sediment inflow to the canal system, however, de- 
pends not only on the sediment load of the stream from which the water is 
diverted, but also on the division of sediment load as the result of flow diver- 
sion. The latter differs considerably according to the local conditions; and 
the data available is so meager that not even a rule of thumb can be made at 
the present. This situation may be improved considerably if data on the 
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division of sediment load under a variety of flow conditions can be collected 
either by taking direct field measurements of sediment load or by applying 
the Einstein bed-load function to both the stream and the canal, after the 
latter has become stable and provided the bar resistance of canal flows has 
been properly described. In many irrigation systems, provisions have been 
made at or close to the canal entrance to eject the coarse material from the 
canal. Today little is known as how effective are these ejectors. In case 

they do exclude all the coarse sediment, the material drawn from the stream 
to the canal will consist mainly of the fine particles. The distribution of fine 
sediment in the stream is generally more uniform, and the division of these 
particles is probably in proportion to the division of discharge at the canal 
entrance. This condition permits the determination of sediment inflow to the 
canal system with ease. Unfortunately there is one more complication. The 
canal inflow no doubt carries with it its share of the wash-load material. 
Since the velocity in the canal is lower than that in the stream, part of the 
wash load may settle out in the canal and form a significant part of the bed. 
Or, in other words, these fine particles which behave like wash-load in the 
stream may behave like bed-material load in the canal, and their transport 
rate must also be included in the stability consideration. How much of this 
part of the load is diverted depends not only on the flow, but also on the 
availability of the material from the upstream watershed. For a small water- 
shed it may be possible to establish a rating curve of the wash load material 
by tedious sediment load sampling program. This task becomes increasingly 
complex if the watershed is large and if the topography, soil, and vegetation 
differ considerably from place to place within the watershed and if the local 
rainfall patterns cover only part of the entire area. On the other hand, it is 
to be remembered that one does not need to know the sediment inflow very 
accurately, and perhaps complications like these can be ignored with little 
consequences, No experience exists on this part, however. 
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Fig. 4 Sediment Rating Curves at Constant Slopes 
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Fig.5 Junction of Stream with a Sizable Tributary 


50,000 cfs 
= 50,000 tons/day 
500 ft. 


@ 


Q2 = 25,000 cfs 
25,000 tons/day 
P2 =(25,000/50,000) ® x 500=350 ft 


Fig. 6 Flow at River Fork 
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